Partially incoherent gap solitons in Bose-Einstein condensates 
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We construct families of incoherent matter-wave solitons in a repulsive degenerate Bose gas 
trapped in an optical lattice (OL), i.e., gap solitons, and investigate their stability at zero and 
finite temperature, using the Hartree-Fock-Bogoliubov equations. The gap solitons are composed 
of a coherent condensate, and normal and anomalous densities of incoherent vapor co-trapped with 
the condensate. Both intragap and intergap solitons are constructed, with chemical potentials of 
the components falling in one or different bandgaps in the OL-induced spectrum. Solitons change 
gradually with temperature. Families of intragap solitons are completely stable (both in direct sim- 
ulations, and in terms of eigenvalues of perturbation modes), while the intergap family may have a 
very small unstable eigenvalue (nevertheless, they feature no instability in direct simulations). Sta- 
ble higher-order (multi- humped) solitons, and bound complexes of fundamental solitons are found 
too. 

PACS numbers: 03.75.Lm; 03.75.Kk; 42.65.Tg; 05.45.Yv 



I. INTRODUCTION 



An ultracold Bose gas is, in general, a mixture of a coherent Bose-Einstein condensate (BEC) and an incoherent 
(fluctuating) "vapor". As shown in Refs. 0,12, the Hartree-Fock-Bogoliubov (HFB) description of such a mixture in 
one dimension is provided by the Gross-Pitaevskii (GP) equation for the condensate order parameter (j){x,t), coupled 
to equations for components of the vapor wave function, u{x,t) and v{x,t), which are responsible for normal and 
anomalous densities of the fluctuations. This approach and related formalisms for the description of fluctuations make 
it possible to analyze various effects, such as quantum phase diffusion in BEC and its depletion in a time-dependent 
trap through transfer of atoms to non-condensed states 's'l, depletion of dark solitons quantum-noise squeezing 
of gap solitons in a repulsive condensate trapped in an optical lattice (OL) deviations from one-dimensionality 
@ , and friction and diffusion of solitons in a cloud of thermal atoms . The full system of coupled HFB equations 
was used in Ref. to show that, in the case of attraction between atoms, the matter flux from the condensate to 
the vapor may lead to splitting of bright matter-wave solitons (which have been created experimentally in ^Li |^ and 
^^Rb 9] condensates) into two fragments that may be regarded as partially incoherent solitons, i.e., bound states of 
the coherent condensate and vapor components, similar to partially incoherent solitons known in nonlinear optics |lClj| . 

The objective of this work is to find solutions for partially incoherent gap solitons (GSs) in an OL potential, and 
investigate their stability, at zero and finite temperatures. GSs in BEC were predicted using the GP equation 
and then created experimentally in a ^^Rb condensate (these solitons were quite "meager", each consisting of 
just a few hundred atoms). Partially incoherent lattice solitons at finite temperature (T) were qualitatively predicted 
in Ref. 'l3l, which relied upon simulations of the GP equation, starting with a random Bose distribution at finite T 
and gradually forming a soliton by switching the OL potential on. Families of two-component GSs of intragap and 
intergap types (see below), i.e., with chemical potentials of the components falling in one or different bandgaps of 
their common OL-induced linear spectrum, were found within the framework of GP equations coupled by repulsion 
between the species 0]. It is relevant to mention that two-component solitons of the intergap type as well as 
their discrete counterparts |l6l |. were earlier predicted in lattice models of nonlinear optics |10|; however, these were 
coherent objects and, unlike the model considered here (see below), they were found in models that did not include 
coherent four-wave-mixing terms. 

The GS is not a ground state of the repulsive condensate (it is obvious that it does not realize an absolute energy 
minimum for the self-repulsive condensate loaded into the OL potential, with a given number of atoms), but it is 
nevertheless stable. It is relevant to stress the difference of the analysis of quantum fluctuations around the GSs from 
the problem of quantum depletion of dark solitons, where fluctuations fill out the notch at the center of the soliton 
and thus gradually destroy it ^4J. For GSs (which are bright solitons), the notch is absent in the family of solutions 
found in the first bandgap [see Fig. [21[a) below]. In the second gap, notch(es) may be present in decaying wings of 
the soliton's waveform [see Figs.|2Ib) and[S], but the soliton's identity is not predicated on them. 

The paper is organized as follows. In Section II, we formulate the model, which is based on time-dependent HFB 
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equations for the coherent condensate and incoherent vapor components interacting with it (at finite T) . Basic results 
for the partially coherent GS families of the intra- and intergap types, found in a numerical form at T = and T > 0, 
are reported in Section III. Examples of higher-order (multi-humped) solitons of various types are presented in Section 
III too. Stability of these solitons is investigated, by means of computation of the corresponding eigenvalues for small 
perturbations, and in direct simulations, in Section IV. Section V concludes the paper. 



II. THE MODEL 



Coupled time-dependent HFB equations are obtained as a truncation of a hierarchy of approximations !l| developed 
for the description of the dynamics of interacting condensate and vapor components of the degenerate bosonic gas, at 
very low but, generally, finite T. Following Ref. 0, the equations for the gas with repulsion between bosonic atoms 
are cast in the following normalized form, with the nonlinearity coefficient scaled to unity: 
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Here L = -(1/2)92^ - ecos(2x) -f 2 (^|0| -I- nj , h = {1 + Afs) \v\ + AfsH'^ and m = - (H- 2Mb)uv* are the normal 
and anomalous fluctuation densities (the asterisk stands for the complex conjugation), 
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is the Bose occupancy, with E the ground-state energy of the tight transverse confinement (note that the one- 
dimensional equations are derived from their 3D counterparts, assuming strong confinement in the transverse plane, 
by means of various approaches implying averaging in the transverse plane [S El), and e is the strength of the 
longitudinal OL potential, whose period is scaled to be tt. A dynamical invariant of Eqs. H1I2|) is the total number of 
atoms, 
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(note that the expression for N explicitly depends on T via Mb)- 

Partially incoherent GSs are looked for as bound states in which the coherent condensate and incoherent vapor 
components are trapped together in the OL, 



(t){x,t) = $(x)e-*^**, {u{x,t),v{x,t)} = {C/(a;) y(a;)e+*'^"*} , 

with chemical potentials and —p,^ subject to the constraint (phase- locking condition), 
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which implies that collisions may kick out pairs of condensate atoms into the vapor (i.e., dynamical depletion of 
the condensate); the above-mentioned transverse-confinement energy was subtracted from the chemical potentials. 
Equations for the stationary parts of the wave functions are obtained by the substitution of expressions (jSJ in Eqs. 
(P) and 101: 
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In the linear approximation, Eqs. l(7|) decouple into the Mathieu equation for the condensate, <&" -I- 
2[fi^ + e cos(2x)] $ = 0, and its rephcas for U and V, giving rise to identical bandgap spectra for the three waves. 
This bandgap spectrum is displayed, for the sake of illustration, in Fig. ^ GS solutions can be found if /x^, Hu 
and jly [the latter chemical potential is equal to — /iu + 2/X0, as per locking condition ^] belong to one or different 
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FIG. 1: The bandgaps (unshaded areas) in the identical spectra of three hn- 

earized equations Q), with fi standing for any chemical potential, fi^, or fly. 

ban dgap s. The respective solutions will be called intragap and intergap solitons, as in Ref. 0| (see also Refs. 
and lig)- In particular, in the typical case of a moderately strong OL, which is used below to present generic results, 
with OL strength e = 5 (i.e., 2.5 Eij, where the recoil energy is, in physical units, Eji = h'^k'^/2m, k being the OL 
wavenumber), the first and second bandgaps in Fig. Qare —2.5 < ^ < 0.7 and 1.2 < ^ < 3.8. 

To conclude the description of the model, it is relevant to stress that the fluctuations, u{x) and v{x), must include 
contributions from all Bloch bands of the periodic potential. For sufficiently large e, the lowest bands are narrow, 
hence the Bloch states in each of them may be approximated by a single mode in u and v. For instance, inspection of 
the spectrum from Fig. Qlfor e — 5 (this value will be used below) demonstrates that the three lowest bands are indeed 
sufficiently narrow to be approximated by single modes, while other bands have little relevance as they correspond 
to very high values of the chemical potential. Furthermore, in this situation contributions from mode-mixing cross 
terms to quadratic quantities (integrated densities, that measure the strength of the fluctuation components in the 
Bose gas, see below) are negligible, in view of the effective mutual incoherence of the Bloch wave functions in distinct 
narrow bands separated by wide gaps. Therefore, in such a representation, the integral quantities actually take the 
familiar form of diagonal sums over several fluctuation modes Q-B- 



III. NUMERICAL RESULTS: FAMILIES OF PARTIALLY INCOHERENT GAP SOLITONS 



A. Intragap solitons 

Generic examples of solutions to Eqs. O with T — and T > 0, in the form of intragap GSs with the three chemical 
potentials falling in the first or second bandgap, are shown in Fig. |21 In terms of Ref. 0|, they are categorized as 
tightly and loosely bound solitons, respectively. Note that A/b = 1 corresponds to T « 70 nK, if the transverse 
trapping frequency is 27r x 1 KHz. In Fig. [3 it is seen that the soliton does not suffer drastic changes with the increase 
of temperature. 

Families of partially incoherent GSs can be characterized by fractions of the vapor components in the total number 
of atoms, i.e., Nu/N and Ny/N [see Eq. Q], as functions of /i„ and fiy, each chemical potential varying within a 
given bandgap. In Fig. |3| we display these dependences for T = 0, with both fiy and jly belonging to the first or 
second bandgap [then /i^ = (/i^ + /i„) /2, see Eq. ©, lies in the same gap, hence the families are of the intragap 
type, indeed]. A "valley" in the plots running along the diagonal means that the symmetric solutions, with fi^ = P-v, 
amount to the ordinary GSs, with u = v = 0. All solitons belonging to these families in the first and second gaps 
feature, respectively, tightly- and loosely-bound shapes, similar to those in Fig. [21 

Figure0]shows that the dependence of Nu/N and Ny/N on temperature is very weak, at fixed values of and jiy 
(in the temperature range considered). Counterparts of the latter dependence for fixed N (rather than fixed chemical 
potentials) may be interesting too, but they need a large pool of numerical data and will be reported elsewhere. 
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FIG. 2: (Color online) Typical examples of stable intragap solitons found for T = and T > 0, which correspond to 
Bose occupancy Mb — and 1, respectively [see Eq. Q] in the first (a) and second (b) bandgaps, with (a) /i„ = —2, 
fly = —2.5, = —2.25; (b) /i„ = 1.5, jly = 3, /i^ = 2.25. In this and other examples, the OL strength is £ = 5. 

B. Intergap and higher-order soUtons 

A family of intergap solitons has been constructed too, with the chemical potentials of the fluctuational components, 
and /lu, belonging to the first and second bandgaps, respectively [then, the condensate's chemical potential, locked 
to /i„ and jljj as per Eq. jnj, /i^ = (/in + fiv) /2, belongs to the second bandgap too]. A typical soliton of this type is 
shown in Fig.[Sl Note that its U{x) component, which corresponds to the chemical potential, which belongs to 
the first bandgap, features a typical tightly bound shape, cf. Fig. E^a), while the shapes of the other two components, 
<^{x) and V{x)^ which pertain to the chemical potentials, /i^ and jl^, which belong to the second bandgap, are weakly 
bound, cf. Fig. Ela). 

The Nu/N and N^/N characteristics for the entire intergap family are presented in Fig. Note that intergap 
solitons with u = v = do not exist, unlike their intragap counterparts, therefore these plots do not feature "valleys", 
unlike Fig. |31 

In addition to the fundamental (single-humped) GSs, various types of higher-order multi-humped states have been 
found too. In Fig. [7| we display (for T — and T > 0) the simplest among them, which is single-humped in the 
condensate field (0), and doubled- humped in one of the vapor components. 

Obviously, the higher-order soliton shown in Fig. His not a bound state of fundamental solitons. On the other 
hand, straightforward bound states can be found too, see an example of a three-soliton complex in Fig. |S1 Note 
that, as per a general principle for the stability of bound solitons on lattices 19], this complex may be stable because 
the phase difference between the bound solitons is tt. As well as the fundamental intra- and intergap solitons, their 
higher-order counterparts of various types form families which fill out the bandgaps. 

IV. STABILITY ANALYSIS 

The stability of the GSs was first tested in direct simulations, which has demonstrated that they are all appear to 
be stable, both at T = and T > 0. Most accurate information about the stability can be obtained from computation 
of eigenvalues for small perturbations, using equations and (|2J) linearized around stationary solitons. In particular, 
the stability of ordinary GSs was previously shown in the framework of the Bogoliubov-de Gennes equations, which 
are derived by the linearization of the GP equation about the solitons 20]. 

Following this approach, we first consider the stability of the GSs with u — v — (i.e., the subfamily along the 
diagonal "valleys" in Fig. against small vapor perturbations, u = e~'^^^ui{x) and v = e^^^'^'*'~^^*vi{x), where A 
is the perturbation eigenvalue. The instability implies the existence of eigenvalues with Im(A) > 0. In this case. 
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FIG. 3: (Color online) The vapor fractions N.^/N and Ny/N in families of stable tightly and loosely 
bound intragap solitons in the first (a,b) and second (c,d) bandgaps versus fiu and fly, at T = 0. 
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FIG. 4: (Color online) The vapor fractions Nu/N and Ny/N vs. temperature [shown 
via Bose occupancy Mb, see Eq. Q] for the same solitons as in Fig. |2] 
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FIG. 5: (Color online) An intcrgap soliton at T ~ 0, with chemical potentials /i„ — 
and /2u = 3, /i^ = 1.5, which fall in the first and second bandgaps, respec- 
tively. The inset illustrates stability of the perturbed intergap soliton in direct simulations. 




FIG. 6: (Color online) The vapor fractions N^/N and N^/N in the family of intergap solitons, at T = 0, versus chemical 
potentials of the two fluctuational components, and /ii,, belonging to the first and second bandgaps, respectively. 
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FIG. 7: Stable higher-order solitons in the first bandgap, for the same parameters as in Fig. OJa). 
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FIG. 8: (Color online) An example of a stable bound state of three fundamental solitons belonging to the first 
bandgap (for /i„ = —2.0, fi^ — —2.2, /i^ — —2.1): the condensate component (a) and vapor wave functions (b). 
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FIG. 9: The instability growth rate for the family of intergap solitons (at T = 0), 
as found from numerical solution of the linearized equations for small perturbations. 



linearized equations (|2l decouple from Eq. yielding 

(l/2)u"+ [ecos(2j:) - 2/9|$ 



Ml + <I>^Wl = Aui, 



l/2)i;"+ ecos(2j;) - 2p|$|^ vi + ^'^ui = -(2^^ + A)wi. 



(8) 



Solving these equations (which do not depend on temperature) numerically (with proper boundary conditions), we 
have concluded that all GSs with zero vapor components are stable against "vaporization" . 

Then, we performed the linear-stability analysis for the full GSs, including nonzero vapor components. We have 
found that the families of intragap solitons in both (first and second) bandgaps are completely stable (for T — Q and 
T > alike), in complete accordance with direct simulations. Preliminary considerations of higher-order intragap 
solitons, such as ones displayed in Figs. [7|and|Hl suggest that they are stable too. 

For the intergap family, a weak instability is revealed by the computation of eigenvalues, see Fig. |51 (since intergap 
solitons cannot exist without vapor components, this instability is specific to the partially incoherent GSs). However, 
this weak instability does not manifest itself in direct simulations (as shown, for instance, by inset in Fig. |SJ), which 
suggests that the intergap solitons, even though being formally unstable, may be observed in experiments. 
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It is also relevant to check the stability of the solitons against deviations from the one-dimensionality [g- In the 
simplest approximation, this amounts to modifying the equations, keeping them effectively one-dimensional but adding 
quintic terms to the cubic nonlinearity 0, ■ Preliminary analysis shows that no additional instability emerges in 
this way. 

V. CONCLUSION 

We have found families of matter-wave gap solitons (GSs) in the degenerate Bose gas with repulsive interactions 
between atoms, trapped at zero or finite temperature in a periodic optical-lattice (OL) potential. Stability of the GSs 
was studied too. The solitons include a coherent condensate wave function, and two components of the incoherent 
"vapor" (which actually comprise many fluctuation modes, due to the OL's band structure). Chemical potentials of all 
constituents of the GS must fall in bandgaps. Accordingly, families of intra- and intergap solitons (including higher- 
order ones, and bound states of fundamental solitons) were found in the two lowest bandgaps, and it was concluded 
that they do not change drastically with the growth of temperature. While the intragap GSs are completely stable, 
their counterparts of the intergap type feature a very small unstable perturbation eigenvalue, but, nevertheless, they 
do not feature any tangible instability in direct simulations. 

We acknowledge valuable discussions with A. Vardi. This work was supported in part by Israel Science Foundation 
(Center-of-Excellence grant No. 8006/03) and U.S. -Israel Binational Science Foundation (grant No. 2002147). 
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